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We present the 12 CO(1-0) and 13 CO(2-l) line maps obtained observing with 
the SEST a ~ 1° x 1° region of the Vela Molecular Ridge, Cloud D. This cloud 
is part of an intermediate-mass star forming region that is relatively close to the 
! Sun. Our observations reveal, over a wide range of spatial scales (from ~ 0.1 to a 

few parsecs), a variety of dense structures such as arcs, filaments and clumps, that 
are in many cases associated with far-IR point-like sources, recognized as young 
stellar objects and embedded star clusters. The velocity field analysis highlights 
i—i ', the presence of possible expanding shells, extending over several parsecs, probably 

q ! related to the star forming activity of the cloud. Furthermore, the analysis of 

the line shapes in the vicinity of the far-IR sources allowed the detection of 13 
molecular outflows. Considering a hierarchical scenario for the gas structure, 
a cloud decomposition was obtained for both spectral lines by means of the 
CLUMPFIND algorithm. The CLUMPFIND output has been discussed critically 
and a method is proposed to reasonably correct the list of the identified clumps. 
.{ph ! We find that the corresponding mass spectrum shows a spectral index a ~ 1.3 4- 

2.0 and the derived clump masses are below the corresponding virial masses. The 
mass-radius and velocity dispersion-radius relationships are also briefly discussed 
for the recovered clump population. 
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1. Introduction 

Vela Molecular Ridge is a complex located in the southern sky, extending in galactic 
longitudes i = 260° 4- 273°, and confined to latitudes b = ± 2°, as reported by May et al. 
(1988), who carried out an extensive survey of the third galactic quadrant (£ = 210° -r- 279°, 
b = ±5°) in the 12 CO(1-0) transition, with a spatial resolution of 0.5°. This complex is 
composed of at least three molecular clouds, named A, C and D (Murphy & May 1991), at 
a distance of d = 700 ± 200 pc (Liseau et al. 1992), and a more distant cloud, named B, at 
d ~ 2000 pc. All these regions (hereinafter, VMR-A, -B, -C, -D, respectively) contain a gas 
mass exceeding a total of 10 5 M Q and are active sites of star formation, as revealed by the 
IRAS mission (Liseau et al. 1992). 

VMR-D is one of the nearest massive star forming regions and therefore one of the most 
studied. A detailed analysis aimed to investigate and describe the star formation activity 
in this cloud has been carried out in the last decade through multifrequency observations 
(Lorenzetti et al. 1993; Massi et al. 1999, 2000; Giannini et al. 2005; Massi et al. 2006a), 
that have pointed out the evidence for young stars clustering. 

At radio wavelengths, further observations of VMR-D have been carried out by Wouter- 
loot & Brand (1999, hereinafter WB99) who mapped, in 12 CO(1-0), 13 CO(1-0), C 18 O(l-0) 
and CS(2-1) transitions, nine zones towards bright infrared sources in VMR field including 
IRS 17, IRS 19, and IRS 20 (Liseau et al. 1992, here we shall use their nomenclature for the 
protostellar IRAS sources in common with their sample). In particular, in the case of the 
three mentioned sources, the association with a gas clump and a molecular outflow has also 
been ascertained. 

Yamaguchi et al. (1999, hereinafter YA99) obtained large-scale 12 CO(1-0) and 13 CO(l- 
0) maps of the whole VMR with NANTEN telescope, with a grid spacing of 8'. They 
also mapped the densest regions in the C 18 O(l-0) line with a grid spacing of 2'. The total 
mass of the complex estimated by means of these three molecular tracers is M gas ( 12 CO)~ 
5.6 x 10 5 M , M 9as ( 13 CO)~ 1.4 x 10 5 M , and M 3fls (C 18 0)~ 5.3 x 10 4 M , respectively. 
Furthermore, these authors also related the gas distribution to the position of the IR sources 
in the field, in order to investigate the nature of the star formation in this complex. 

As a further progress in the knowledge of the VMR gas morphology, Moriguchi et al. 
(2001) repeated the 12 CO (1-0) observations of the ridge using the same instrumentation of 
YA99, but mapping it with a finer grid spacing of 2'. These authors remark the highly 
filamentary distribution of the gas emission and discuss the possible interaction between 
VMR and the well-known Vela supernova remnant (SNR), associated with a pulsar located 
at t ~ 263.5°, b ~ —2.5°. They conclude that the molecular gas is associated to the SNR, 
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but it is pre-existent. 

More recently, Fontani et al. (2005), studying a conspicuous sample of high- mass pro- 
tostellar candidates in the southern sky, considered also IRS 21 and found that this object 
is associated with CS(2-1), CS(3-2), C 17 O(l-0) lines and 1.2 mm continuum emission peaks. 

Dust continuum millimeter observations ofa~l°xl° portion of VMR-D (the same as 
investigated in this paper) have been obtained with SIMBA/SEST (Massi et al. 2005; De 
Luca et al. 2005; Massi et al. 2006b) revealing both filamentary and clumpy appearance. In 
particular, De Luca et al. (2005) compiled a list of 29 detected clumps (discussed in Massi et 
al. 2006b, hereinafter MA06) in most cases exhibiting non-gaussian shapes and a tendency 
to clusterize, both suggesting a multiple star formation activity. 

In addition to this millimeter survey, further SIMBA observations are available for the 
neighborhoods of IRS 17 (Faiindez et al. 2004; Giannini et al. 2005), IRS 16, IRS 18, and 
IRS 20 (Beltran et al. 2006), that allowed to detect dust cores and determine their masses 
and luminosities. 

The present work contributes to investigate the star formation activity in a significant 
part of VMR-D by means of molecular line observations down to spatial scales smaller than 
those investigated in the preceding literature. This has been possible by exploiting the spatial 
resolution of the SEST telescope that has been used to observe this cloud in two millimetric 
bands corresponding to the transitions 12 CO(1-0) and 13 CO(2-l), the first one suitable for 
tracing the large-scale distribution of the gas, the second one useful for better probing the 
highest column densities. In this way, several aspects related to the star formation can be 
investigated, such as, e.g., the gas distribution and kinematics to scales smaller than 1 pc (the 
typical size of pre-stellar cores), the correlation between the intensity peaks of the maps and 
the position of the brightest infrared point sources, the presence and census of proto-stellar 
outflows in this region. 

In this paper we focus our attention on the morphological and kinematical aspects of the 
gas distribution. In Section 2 we provide a short description of our observations and of the 
data reduction procedure. In Section 3, after presenting the integrated intensity and channel 
maps of the emission in both lines, we discuss some significant velocity-position (vel-pos 
hereinafter) diagrams, which highlight the complex structure of the gas velocity field in this 
cloud. Section 4 illustrates the search for clumps in the cloud, based on the CLUMPFIND 
code, and presents the derived mass distribution. Finally, in Section 5 a summary of the 
main results is given. 
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2. Observations and data reduction 

In the cloud identification by Murphy & May (1991), VMR-D is located in the rightmost 
end of the galactic longitude interval of the complex, and shows two main CO emission peaks. 
The same characteristic is recognizable, but in greater detail, in the 12 CO(1-0) and 13 CO(l- 
0) maps of YA99 and in the 12 CO(1-0) map of Moriguchi et al. (2001). Our observations 
encompass the region of the higher longitude peak (the smaller in size), located at I ~ 
263° + 265°, b ~ + -1° (see also Figure 1 in Massi et al. 2003). Hereinafter, with VMR-D 
we shall indicate this part of cloud. 

VMR-D was mapped in the 12 CO(1-0) [y = 115.271 GHz) and 13 CO(2-l) {y = 220.399 GHz) 
lines with the 15-m Swedish-ESO telescope (SEST, see Booth et al. 1989) at La Silla, Chile, 
during two complementary observational campaigns, in September 1999 and January 2003, 
respectively. 

A high-resolution Acousto-Optical Spectrometer was used as a backend, with a total 
bandwidth of about 100 MHz and a resolution of 41.7 kHz, in splitted mode, i.e. in a 
configuration allowing the use of both 115 and 230 GHz receivers simultaneously (1000 
channels for each receiver). This corresponds to velocity resolutions of ~ 0.12 km s _1 at 
115 GHz and ~ 0.06 km s _1 at 220 GHz. The FWHM of the primary beam and the main 
beam efficiency factors are 45" and 0.7 at 115 GHz, and 23" and 0.5 at 220 GHz, respectively. 
During the 1999 run only, the 230 GHz receiver was tuned to the frequency of 13 CO(2-l) and 
the corresponding mapped area is then smaller than in the 12 CO(1-0) line. 

All the observations have been obtained in the frequency switching mode (Liszt 1997), 
with the switch interval depending on the line frequency. This has been chosen to be larger 
than the extent of the emission (to avoid line signal overlap), but small enough to avoid 
losses of emission features in the reference cycle, even if this was not always the case for the 
13 CO(2-l) observations, in which the spectral range is halved with respect to the 12 CO(1-0). 

We acquired the maps adopting a grid spacing of 50", so that the 12 CO(1-0) is slightly 
undersampled while the 13 CO(2-l) is a factor of two undersampled. For an estimated distance 
d = 700 pc, the adopted grid spacing corresponds to a spatial scale of 0.17 pc on the 
cloud. The coordinates of the (0,0) position in the map are a(2000) = 8 ft 48 m 40 s , 5(2000) = 
— 43°46'12", and the pointings range from a jf = —2000" to a fj = +1200" in right ascension 
and from 5 jf — 2000" to 5 jf = +1950" in declination. Some spatial gaps, mainly in the 
case of the 13 CO(2-l) line, remained after our observing campaign and are described below. 

The integration time at each point was generally set to t int = 10 s, but there is a 
significant fraction of spectra observed with t int = 20 s. The typical rms noise affecting the 
data (main beam temperature) is AT rms ~ 0.7 K for both lines. 
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Towards many pointings, the observations were repeated on different dates in order to 
check the data for consistency. In particular, the matching between 1999 and 2003 spectra 
was fully satisfactory (differences between peak temperatures are always within 10%), except 
for a raster of points located at south-east of the map, in the +250" < a Q ff < +950", 
—2000" < S ff < —1050" region, which were affected by temporary instrumental problems, 
and then partially re-observed on subsequent dates. The reader should then be aware that 
in the region +250" < a off < +950", -2000" < S off < -1700", these repeated pointings 
are not available, so this part of the map is less accurate. 

The pointing accuracy was checked every 2-3 hours towards nearby (in the sky) SiO 
masers being within ~ 5". 

Data reduction followed the pipeline described in Massi et al. (1997): first, spectra 
in antenna temperature were scaled by the main beam efficiency factor 7]^, in order to 
express them in terms of main beam temperature (T m t> = TA/i] m b); then a polynomial fit of 
the baseline was subtracted and a folding was performed on the resulting spectra. Finally, 
spectra have been resampled to a velocity resolution of 0.12 km s _1 for the 12 CO(1-0), and 
0.06 km s- 1 for the 13 CO(2-l), respectively. 

In the cases of spatial superposition of repeated observations, the corresponding spectra 
have been averaged (with weights depending on the integration time and the inverse of 
the system temperature), in order to obtain a better signal-to-noise ratio for the resulting 
spectrum. 



3. Observational results 

3.1. Integrated intensity maps 

Out of the 4258 observed points, significant 12 CO(1-0) emission was detected from 3392 
of them, corresponding to a detection rate of ~ 80%. In Figure 1, a sample of reduced 
12 CO(1-0) spectra, taken towards the north-western region of our VMR-D map, is shown. 
Despite this part of the grid is relatively small, it is representative of the velocity field 
complexity that characterizes this cloud. In the case of the 13 CO(2-l) line the points showing 
significant emission at our sensitivity level are 648 out of 2393 (~ 27%). 

The bulk of the emission falls in the range V/ sr = —2 + 20 km s" 1 for the 12 CO(l- 
0) line and Vi sr = + 14 km s _1 for the 13 CO(2-l) line, as shown in Figure 2, consistently 
with WB99 and YA99. The total intensity of the 12 CO(1-0) and 13 CO(2-l) emission, i.e. 
J T m b dv integrated from -2 to 20 km s _1 , are represented in grey-scale in Figure 3 and 
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Figure 4, respectively. 

The 12 CO(1-0) map in Figure 3 evidences a complex clumpy and filamentary structure 
displaying cavities of various sizes. The analysis of the velocity components, based on the 
channel maps and the vel-pos diagrams presented in Section 3.3, helps us in separating the 
various emitting regions. Here we limit ourselves to sketch the main characteristics of the 
cloud shape, using the locations of the bright far-IR sources in the field as reference points 
(see Figure 5): 

• A large region in the northern part (-1500" < a off < +1200", +1100" < 5 off < 
+ 1950") characterized by a strong, distributed emission. The coincidence between the 
position of IRS 20 and the most intense peak in this part of the map is evident. 

• A central NE-SW region (-2000" < a off < +500", -750" < 6 off < +1000") ap- 
proximately elongated from the location of IRS 19, associated with a bright peak of 
integrated emission, to that of IRS 17 (the strongest peak in the whole map) and 
IRS 16, at the western boundary. This emitting region separate three zones, east, 
west, and south of it, respectively, characterized by an evident lack of emission and 
determining an arc-like appearance for the gas. 

• An apparently compact structure in the south-eastern part of the map (—250" < 
a off < +1050", -1800" < S off < -800"), hosting IRS 21. This region contributes also 
to form the southern arc of molecular gas. 

The distribution of 13 CO(2-l) emission closely follows that of the more intense 12 CO(l- 
0), as can be easily seen in Figure 4 and especially in Figure 5, in which, to highlight the 
correlation between the emission of the two lines, we superimpose the levels of the latter on 
the grey-scale map of the former. The 13 CO(2-l) traces the densest parts of the 12 CO(l- 
0) map, and its peaks coincide with those of the 12 CO(1-0), except for the main peak in 
the south-eastern region. The northern region of the map is characterized by a number of 
isolated peaks, some of which are also recognizable in the 12 CO(1-0) map. 

In Figure 5 the positions of 25 IRAS PSC sources with fluxes increasing with wavelength 
(F12 < -F25 < Fqo) are also overplotted. Among these objects, there are also IRS 16, IRS 17, 
IRS 19, IRS 20, IRS 21, characterized by the further constraint F 2 $ > 2.5 Jy, recognized as 
intermediate mass YSOs by Liseau et al. (1992), and associated to embedded young clusters 
(Massi et al. 2000, 2003). Incidentally, a search for objects with increasing fluxes in the 
MSX Point Source Catalog (Egan et al. 2003) provided only four sources, coinciding with 
IRS 16, IRS 17, IRS 19, IRS 20. The reason for this small number of selected objects with 
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respect to the IRAS archive query probably lies in the different extension of the spectral 
range observed. 

Typically, all these sources are found in the densest parts of the cloud, and in most 
cases their locations follow the arc-like morphology of the gas. This is suggestive of a star 
formation process triggered by the effects of expanding shells driven, e.g., by nearby massive 
young stars or supernova remnants. In Section 3.4 we shall explore in some detail this 
possibility. 



The present observations have been also used to derive some physical properties of the 
emitting gas. As a first approximation, the column density of the molecular gas traced by the 
thermalized and optically thick 12 CO(1-0) line can be determined by means of the empirical 
formula 



which is based on a galactic average factor (Strong et al. 1988). The molecular mass can be 
derived by the relation 



(see, e.g., Bourke et al. 1997), where d = 700 pc is the distance of VMR-D, (d 2 AaA5) is the 
size of the emitting area at the observed position, ran is the mass of a hydrogen atom, and \i 
represents the mean molecular mass. Adopting a relative helium abundance of 25% in mass, 
fj, = 2.8 and the derived cloud mass is M 9as ( 12 CO) = 1.5 x 10 4 M . The quite large relative 
error on the distance estimate, Ad = 200 pc, dominates the total uncertainty affecting this 
value (AM/M ~ 50%), although it is necessary to clarify that the errors implied in the 
assumptions underlying this calculation probably exceed and dominate this amount. This 
value represents only a small fraction (~ 2.5%) of the mass reported by YA99 for the whole 
complex, i.e. an area of ~ 100 deg 2 , because this is the mass of only a small part of one 
out of four clouds. Note also that VMR-B strongly contributes to the determination of the 
mass of the whole complex (see Equation 2), being three times more distant than the other 
clouds and showing a similar angular extension. 

A different and more rigorous method to derive the column density relies on the simulta- 
neous observations of the two lines, exploiting their different optical depth. As a first step the 
excitation temperature is estimated from the peak temperature of the 12 CO(1-0) emission, 



3.2. Gas physical properties 




(1) 




(2) 
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indicated with T R ( 12 CO), assuming LTE conditions and r » 1 (e.g., YA99): 

5 ' 53 

ex ~ ln{l + 5.53/[T R ( 12 CO) + 0.819]} ' ^ ' 

The peak strengths were derived by gaussian fits to the 12 CO(1-0) line profiles. Often, these 
show very complex profiles in this region: in these cases, we fitted only the component corre- 
sponding in velocity to the 13 CO(2-l) feature, that generally shows only a single component. 
Assuming that the excitation temperature is the same also for 13 CO we then determined the 
optical depth of the 13 CO (2-1) per velocity channel by means of the relation: 

r ^ ( i3 C0 ) = _ In { 1 - T - (13CQ) ( ^ 0.02^1 ~* 1 (4) 

vK } \ 10.58 Vexp(10.58/T e:c ) - 1 J J 1 ; 

and the corresponding column density by means of the relationship (e.g., Bourke et al. 1997) 

M(^m\ i 91 v i ni4 ( T ex + 0-88) exp(5.29/T ea ) f 13 

N( CO) = 1.21 x 10 - _ exp( _ 10 58/Tex) J r v ( CO) dv (5) 

Finally, after N(H 2 ) is derived from N( 13 CO) assuming a CO abundance ratio of 
7 x 10 5 (Dickman 1978), the molecular mass is calculated in the same way as for 12 CO(1-0), 
amounting to M 9as ( 13 CO) = 1.2 x 10 3 M . 

Given the undersampling of 13 CO(2-l), a twofold bias (acting in different directions) 
affect the mass determination obtained through this line. On the one hand, some emission in 
this line is undoubtedly lost and hence unaccounted for in the total mass determination, but 
on the other hand, being 13 CO(2-l) less beam-diluted than 12 CO(1-0), then the excitation 
temperature as derived from 12 CO(1-0) peak main beam temperature is underestimated 
(because of the smaller beam filling factor), causing an overestimate in the derived column 
densities. 

In order to compare the LTE method results with those obtained by applying Equa- 
tion 2, we calculated the mass traced by 12 CO(1-0) considering only those positions showing 
detectable 13 CO(2-l) emission; we found M; as ( 12 CO) = 4.8 x 10 3 M . 

Both the M ffas ( 12 CO) and M 9as ( 13 CO) masses are comparable to those obtained for the 
Orion molecular cloud which possesses similar characteristics from the point of view of the 
star formation activity. Of course, the comparison has to be performed taking into account 
both different angular size and distance of the two clouds. Mapping in 12 CO (1-0) an area 
of 29 deg 2 in Orion A and of 19 deg 2 in Orion B, both located at 560 pc (Genzel & Stutzki 
1989), Maddalena et al. (1986) derived M A = 1 x 10 5 M and M B = 9 x 10 4 M , respectively. 
Subsequently Cambresy (1999) derived, from a visual extinction analysis of the whole cloud, 
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a total mass M AB = 3 x 10 5 M . Finally, a mass estimation of Orion A through 13 CO 
emission, although considering the J = 1 — transition, has been provided by Nagahama 
et al. (1998), that found M Ajl3 = 5.4 x 10 4 M Q . Considering then the observed areas, these 
values appear consistent with those estimated for VMR-D. 

3.3. Velocity structure 

In Figures 6 and 7 we present the maps, for both lines, of the integrated intensity taken 
in velocity intervals of 2 km s" 1 and in steps of the same amount. These velocity channel 
maps clearly show how the emission detected in different locations is contributed by multiple 
velocity components. In Table 1 the characteristic velocity ranges for the main areas of the 
maps are summarized. 

The 13 CO(2-l) channel maps, in particular, allow us to better recognize the main velocity 
components of the densest regions. For example, the emission from the northern part of the 
map appears clearly separated into two components: one in the Vi sr ~ + 4 km s _1 range, 
and the other in the V\ sr ~ 6 + 12 km s" 1 range. 

Due to the large number of observed positions, as well as the geometrical regularity of 
the grid (in particular for the 12 CO(1-0) transition), we can extract vel-pos diagrams, once 
an offset in a or 5 is given and the observed spectra along the given strip of the map are 
considered. Here we present 24 declination vs velocity 12 CO(1-0) stripes with a fixed a ff 
spanning the range from a ff = +950" to a Q ff = —1450" (i.e. the location of IRS 17) in 
steps of 100" (Figure 8), that visualize the velocity field across the cloud. In principle the 
same could be done for the 13 CO(2-l) line, but at our sensitivity level the vel-pos diagrams 
do not add further important information so that here we omit them. 

Many clump-like structures are visible in these diagrams, confirming that VMR-D pos- 
sesses a high degree of inhomogeneity as will be discussed in Section 4. Several cases of strong 
line broadening are also detectable in correspondence to the brightest emission peaks, as for 
example at a Q ff = +50" or a Q ff = —1350", confirming the probable presence of outflows. 

The appearance of arc-shaped structures in these diagrams can be explained either by 
the presence of two Galactic velocity components corresponding to gas located at different 
distances or by the action of expanding shells. In particular, this kind of structure is de- 
tectable in the regions +750" < a off < -150" and -1800" < 5 off < -500" (southern part 
of the map), and +350" < a off < -750" and +400" < 5 off < +1800" (northern part of the 
map). 
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The first hypothesis implies extremely different kinematical distances (~ 1 kpc) between 
the two emitting components, an occurrence indicating the presence of uncorrelated clouds 
along the same line of sight. The other possibility of the shell expansion is intuitively 
suggested by the arc-like shapes (in Figure 8 we draw ellipses to show four possible cases), 
even if this can be also justified by simply invoking the internal motions of the cloud, a 
scenario that seems more appropriate as will be discussed in next section. 



The presence of arc-like features in both a — 5 (Figures 3, 4, 5) and vel-S maps (Fig- 
ure 8) of VMR-D region could be interpreted as a signature of possible expanding shells, a 
phenomenon that is often related with star forming activity. 

To explore this possibility we considered the available information on the distribution 
of the known OB stars and Hn regions in the VMR-D field, that are shown in Figure 9 along 
with the spatial distribution of the gas. Three arc-like molecular structures are marked with 
dashed lines, but only in the case of the eastern arc we could guess an association with clear 
driving sources candidates. A possible projective correspondence can be established between 
the Gum 18 Hn region (Gum 1955) and this arc, that could be tentatively related to it, but 
not for RCW 35 (Rogers et al. 1960). However, the centroids of the two structures do not 
coincide and the distance of Gum 18 is unknown. In addition, the OB stars V* OS Vel and 
CD-43 4690 are not far from the apparent centroid of the molecular arc, but their distances 
(1.7 ±0.3 kpc in both cases, Russeil 2003) are inconsistent with that estimated for VMR-D, 
so that we conclude that they cannot be responsible of the arc-like morphology of the gas. 
Similar considerations apply in the case of the possible association between the southern arc 
and the location of HD 75211 (d = 1.2 kpc, Savage et al. 1985). 

By simple models, describing the expansion of a Hn region and the fragmentation of the 
shocked dense layer surrounding it, it is possible to estimate the timescales involved in the 
hypothesis that the shell driving sources are Hn regions. The relation between the radius R 
of the shell, expanding in a homogeneous and infinite medium, and its lifetime is 



(Spitzer 1978), where Ro is the Stromgren sphere radius (in pc), and c\\ the isothermal 
sound speed in the ionized region (in km s _1 ). Being the molecular arc ~ 1000" in radius, 
the corresponding spatial scale is R ~ 3.4 ±1.0 pc at the estimated distance of VMR-D, 
where the uncertainty is related to the distance. For a star emitting 10 49 ionizing photons 



3.4. Expanding shells 




Myr 



(6) 
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per second, we find Rq = 0.677 n 3 pc, where n.3 = no/(1000 cm -3 ) and no is the ambient 
density. Assuming the sound speed c\\ = 10 km s^ 1 and varying the density in the range 
n = 300 4- 3000 cm -3 , we find t dyn ~ 0.2 4- 1.1 Myr. The scaling of the age with the density 
is shown in Figure 10 for different estimates of the distance. 

Considering now the Whitworth et al. (1994) model describing the fragmentation of the 
molecular gas under the action of an expanding Hn shell, it is also possible to determine the 
time elapsed before the onset of this phenomenon. For the considered star this model gives 

tf ra g = 1.56 c 7 / 11 0/11 Myr (7) 

where c.2 = c s /(0.2 km s -1 ), and c s is the sound speed in the molecular gas, an important 
parameter of the model we chose in the range c s = 0.2 4 0.6 km s -1 . 

Because both fragmentation and star formation appear to be already in progress, the 
relation between these two times have to be td yn > tfrag, and this happens when td yn ^ 
8 x 10 5 yr and no > 1600 cm -3 , if we consider the preferred values d = 700 pc and c s = 
0.2 km s -1 . The effect of a different distance is also shown in Figure 10. 

We note that these timescales are consistent with the age of the cluster formation activity 
in this region, estimated to be within 1 -=-6 Myr by Massi et al. (2006a). This makes possible 
that both arcs and the star formation activity in them are caused by the expansion of Hn 
regions. However, since no clear excitation source candidates have been found in the vicinity 
of the arc centers, we cannot exclude different possible origins for the shells, such as for 
example stellar winds from a previous generation of young stars. 

Despite another possibility is suggested by the well-known Vela SNR, apparently in- 
teracting with the whole VMR on larger scales (see Figure 1 in Moriguchi et al. 2001), a 
correlation with the southern arc in VMR-D is far from clear. Furthermore, considering 
the possible connection between this molecular arc and its star formation activity, there is 
a clear inconsistency between the age of the YSOs and that estimated for the Vela SNR 
(t SNR ~ 10 4 yr, see, e.g., Moriguchi et al. 2001). Two further remnants, SNR 266.3-01.2 
and SNR 260.4-03.4, are too far for being responsible for the peculiar morphology of the 
investigated region. 

In this respect, the hypothesis that the fluidodynamic evolution of the gas can be re- 
sponsible of the observed filamentary condensations can not be rejected. 

As mentioned above, also the analysis of the cloud velocity field can provide useful 
information about the presence and the physical properties of possible expanding shells. 
For example, the shell-shaped structure highlighted in Figure 8 at a a ff = +350" can be 
associated, due to its spatial location, to the southern arc identified in the a — 5 map. It 
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is characterized by a peak separation AVJ sr ~ 7 km s -1 , corresponding to an expansion 
velocity V exp ~ 3.5 km s _1 . Adopting a spherical symmetry, we estimate a mean radius 
i? sh ~ 5.5 pc and a corresponding dynamical age td yn — 1-5 Myr, a value consistent with the 
possible timescales calculated above for an expanding Hn region. A rough estimate of the 
kinetic energy can be also obtained by adopting the assumption inspiring the Equation 2 
and considering the expansion velocity of each spectral component. In this way we find 
E kin ~ 2.3 x 10 47 erg, a value comparable, for example, with the kinetic energy of a Hn 
region possessing the same radius and age (Spitzer 1978). The issue of the identification 
of a driving source for this shell remains open, as already discussed when we examined the 
corresponding arc observed in the a — 5 map. 

Another shell-like structure is well recognizable in the northern part of the vel-pos 
diagrams, best visible at a Q ff = —450". With considerations similar to the previous case, 
we find V^p — 4.5 km s -1 , R s h ^ 7 pc, td yn — 1-5 Myr, and E kin ~ 5.4 x 10 47 erg. Also in 
this case the identification of a possible driving source is uncertain: some IRAS red sources, 
namely those identified as 9, 11, and 14 in Table 2, are located in the neighborhoods of the 
shell, but these objects might be presumably the product of a star forming activity triggered 
by an expanding shell, rather than the cause of it. In conclusion, the ages derived for these 
two candidate shells are similar, and support the hypothesis of a synchronous intermediate- 
mass star formation in this cloud, originated by a strong compression of the gas due to the 
associated shocks. For this process we suggest an upper limit in timescale of t ~ 1.5 Myr. 

In Figure 8 two further arc-shaped structures have been marked with solid lines, around 
a D ff = —650" and a D ff = —750", respectively. Although less evident than the cases consid- 
ered above, these shapes support the hypothesis of a coupling between the arcs recognized in 
the integrated intensity maps (in these cases, the western and the southern one, respectively) 
and those present in the vel-pos diagrams. 

3.5. Search for outflows 

Since several positions in the map show a clear broadening of the line profile (see, e.g., 
the diagram at a = —1350" of Figure 8) we studied these positions in more detail because 
this is a typical signature of the presence of energetic mass outflows (see, e.g., Bachiller 
1996). 

To this aim we considered the locations of the 25 IRAS sources shown in Figure 5 and 
quoted in Table 2, carrying out a systematic search for outflows in the 12 CO(1-0) line toward 
these objects. The detection method is the same as adopted by WB99: a gaussian is fitted to 
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the upper part of a spectral feature, centered on the line peak, and at the same intensity; the 
smaller of the two intervals between the peak velocity and the velocities at half maximum is 
taken as the HWHM of the fitting gaussian. When the profile wings significantly exceed the 
gaussian wings, the presence of an outflow is inferred and the blue and red components are 
obtained by integrating the difference T m b — Tfu between the half maximum of the gaussian 
and the velocity where the line wing fades into the noise. In Figure 11, panels a and d, 
we show, as an example, the 12 CO (1-0) spectrum of the pointings closest to the IRS 17 and 
IRAS 08474-4325 (identified with 3 and 17 in our list, respectively) locations, with the fitting 
gaussian determined as described above. In panels b-c, e-f the intensities of the blue and 
red components, in the neighborhoods of these two sources are represented in greyscale. 

In seven cases the presence of one or more additional spectral components blended with 
the main feature made impracticable a clean detection of any outflow that might be present. 
This difficulty restricts our sample, even if it remains still meaningful in the perspective of 
providing good target candidates for future APEX and ALMA observations. 

In 13 out of the remaining 18 cases, we detected a possible outflow, as summarized in 
Table 2: column 9 contains a flag indicating the result of the detection, column 11 the width 
of the velocity range subtended by the wings, and column 12 the offset of the wing maximum 
emission from the map pixel in which the source falls. Note that in two cases only one wing 
is detected and that when the detection flag "Y" is followed by an asterisk, a gaussian fit of 
a secondary (but well-separated) component has been subtracted to correctly estimate the 
wing intensities. 

Unfortunately, the complexity of the 12 CO(1-0) spectra in this region makes sometime 
difficult an accurate detection of the wings by simply applying the WB99 method. The 
example of IRS 17 in Figure 11, panel a, is illustrative: the temperature peak does not stand 
in the center of the feature, so that a much larger amount of emission is recognized in the 
red wing. Instead, the blue wing seems better suited to be interpreted as an outflow. Vice 
versa in the case of IRAS 08474-4325: the two cases discussed have been shown in Figure 11 
just to highlight the differences between a reliable detection of outflows (blue wing of IRS 17 
and red wing of IRAS 08474-4325), and more controversial cases (red wing of IRS 17 and 
blue wing of IRAS 08474-4325). The latter are indicated in Table 2 with a dagger mark. 

Because our work was not focused on the study of the outflows, we remark that the 
grid spacing of 50" is quite coarse for a detailed description of the outflow morphology, and 
a larger integration time would have been required for each spectrum to better disentangle 
the wing tails from the noise. These problems are evident in comparing our results with 
those obtained by WB99 for the three objects in common with our sample (associated with 
IRS 17, IRS 19, IRS 20 respectively): in our case the wing width is significantly smaller than 



-14- 



found by WB99, ranging from 32% (blue wing of IRS 20) of their values to 71% (red wing 
of IRS 17). 

Furthermore, considering that 13 CO(2-l) observations are too noisy for a similar anal- 
ysis, it is almost impossible to exploit them for a quantitative study of the outflow physical 
properties. In any case our analysis significantly improves the outflow statistics in the ob- 
served region, from three to 13 objects, providing a larger sample of targets for future higher 
resolution observations. 

4. Cloud decomposition 

To characterize the structure of VMR-D and the distribution of its velocity components 
in a quantitative way we applied cloud decomposition techniques to our data, a tool typically 
adopted to study possible hierarchies in the molecular cloud structure. In particular, we 
adopted the 3D version of the CLUMPFIND (hereinafter CF) algorithm (Williams et al. 
1994), running on a a-8-v cube containing all the reduced spectra for each line. Important 
adjustable parameters in the CF scanning procedure are the radiation temperature threshold 
T m i n and the level increment AT, the first one being generally set slightly above the noise 
fluctuations. This is an important point that requires accurate tests because larger and 
smaller thresholds involve the loss of faint peaks and the inclusion of noise peaks, respectively. 
The second parameter AT controls the separation of close peaks, as clearly discussed in Brunt 
et al. (2003). 

A large number of tests have been run on the 13 CO(2-l) data cube, which is characterized 
by a smaller number of features and a reduced diffuse emission, two characteristics that 
allow a better control and verification of the reliability of the clump assignments. To obtain 
a convenient reduction of the noise effect and to ensure some stability to the results, the 
spectra were preliminarily resampled in velocity by a factor of 5. 

Exploring the input parameter space, we noted that the choice of the parameter set can 
produce not only systematic (as those described above), but also random effects. These are 
particularly evident, for example, in large clumps corresponding to a single bright line, that 
are sometimes decomposed in two or more clumps, an effect that is illustrated in Figure 12. 
We verified that, adopting a different set of parameters, the clump could appear as single, the 
price being that some other cases, previously recognized as single, now appear as separated. 
This experience suggests that the tuning of the parameter values is critical with respect to 
the number of detections. This kind of problems in cloud decomposition with CF have been 
extensively discussed in Brunt et al. (2003), Rosolowsky & Blitz (2005), and Rosolowsky & 
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Leroy (2006). 

In this work, the relatively low number of detected clumps suggested us a very simple 
approach based on the careful inspection of a large number of CF outputs. First of all, 
we ran CF adopting a parameter set ensuring a good detection of the features, i.e. a T min 
value as low as possible, although above the noise limit, and a similar value also for AT, to 
minimize the effect of noise fluctuations affecting the spectra. Both values have been set to 
2 K, corresponding to ~ 3cr in our map. In this way, the cataloged emission accounts for the 
66% of the total emission in the 13 CO(2-l) map, as shown in Figure 2. 

The obtained clump list has then been analyzed looking for those objects that, despite 
they are quite close in the a-S-v space, clearly appear as spuriously splitted by CF. We 
merged together these clumps to obtain a final clump distribution as free as possible from 
the most apparent misinterpretations of CF. In this way, we do not claim to determine the 
"true" clump distribution in the cloud, but we believe that the census resulting from our 
clump merging method is more reliable than the simple CF output. 

Our experience in analyzing the clump merging suggested a set of constraints that have 
been included in an automatic merging algorithm. In fact, we merged only clump pairs 
whose centroids satisfy all the following criteria: i) both distances in a and 5 are less or 
equal to 3 grid elements (corresponding to a spatial scale of ~ 0.5 pc); ii) the difference 
between the central velocities is less than or equal to 3 channels (i.e. Av < 0.9 km s _1 , after 
resampling); iii) the a-S-v space total distance d = \J Act 2 + A5 2 + Av 2 is less than or equal 
to 5 volume pixels. In this context the pairs showing a spatial discontinuity in their emission 
have been taken as separated. 

After this procedure the number of detected 13 CO(2-l) clumps decreased from 67 to 
49 clearly modifying the resulting mass spectrum, with the total detected mass remaining 
obviously unchanged. 

The list of these clumps, along with their coordinates, velocity centroids, radii, masses, 
corresponding virial masses and optical depths is presented in Table 3. 

The clump masses have been evaluated assuming the LTE condition, as we did in Sec- 
tion 3.2 for the whole cloud, so that the total mass assigned to these clumps amounts to 
845 M , corresponding to the 77% of the whole M 9as ( 13 CO) calculated. 

The virial masses in column 8 are obtained from the observed radius R and velocity 
dispersion AV under the hypothesis of constant density (MacLaren et al. 1988) 
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For density radial profiles as p oc r 1 and p oc r 2 the multiplicative constant decreases to 
190 and 126, respectively. 

The optical depth is calculated along the line of sight of the clump centroid using 
Equation 4. We list in column 9 the optical depth at peak velocity, and in column 10 the 
integral of this observable over the velocity channels assigned to the considered clump. Note 
that for the clump VMRD43 the maximum and the integrated optical depth are significantly 
higher than others, because of the similarity of 12 CO(1-0) and 13 CO(2-l) peak temperatures. 

For clumps at the boundary of the observed zone a flag "X" is added to indicate that 
the quoted clump mass value actually is a lower limit. 

Comparing these results with the list of 29 dust cores presented in MA06 and obtained 
applying the 2D version of CF to 1.2 mm continuum observations of VMR-D at the SEST, 
we searched for possible associations between 13 CO(2-l) and dust condensations, adopting 
as a criterium the correspondence between the dust core centroid and one of the map points 
assigned by CF to a gas clump. We found 16 cases that are reported in Table 3. Note 
that in this way the information on the velocity field, allowing to separate possible multiple 
components along the line of sight, can induce the association between two gas clumps and 
a group of dust cores (see, for example, the clumps VMRD36 and VMRD37). On the other 
hand, by comparing the two lists in the opposite way, we can note that only the dust cores 
MMS10 and MMS11 are not associated to 13 CO(2-l) clumps. This implies that, in general, 
the continuum emission is able to track even denser zones than 13 CO(2-l) line, and in fact, 
considering the gas clumps associated to a group of MA06 dust cores, we find that the total 
mass is larger or quite similar to the mass traced by dust. The average value of their ratio 
is in fact M gas (line) /M gas (continuum) = 1.44 ± 0.11, while the median of the distribution is 
0.96. Furthermore, the clumps that remain unassociated have intermediate or small masses 
(M < 20 M ). Finally, in Table 3 the associations with the IRAS PSC sources of Table 2 are 
also indicated: 22 out of 25 objects fall in the area observed in 13 CO(2-l), and the locations 
of 15 sources (among which are the five IRS objects of Liseau et al. 1992) can be associated 
with the recognized clumps. 

The derived mass spectrum is shown in Figure 13, with Poisson error bars, along with 
the spectrum resulting from the original CF output. For nine clumps we can only give a lower 
limit to the mass because they are bounded by our map limits. This effect is particularly 
important in the case of the large clump hosting IRS 16 (clump VMRD1 in Table 3). The 
minimum mass detectable and the completeness limits, calculated as in Simon et al. (2001) 
by setting a 10<r confidence level (see also Bains et al. 2006), and estimated in 0.19 M and 
0.98 M , respectively, are also shown. 
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A linear least-squares fit of the distribution of the masses above the completeness limit 
provides an estimate of the dN/dM oc M~ a law exponent, which has been found to be 
a = 2.0 ± 0.3 for the "corrected" sample and a — 1.8 ± 0.3 for the original CF output (see 
Figure 13). In general, dealing with a relatively small sample of clumps, the spectral index 
can be also sensitive to the particular choice of the bin size. It was set here to be 0.2 (in 
logarithmic units), after verifying the stability of the resulting slope with respect to other 
possible choices. By varying, as in Bains et al. (2006), the bin size in the range of ±10% we 
obtain slopes that always fall within the associated error. 

Concluding this discussion about the 13 CO(2-l) clump mass distribution, we remark 
that the observed values are significantly below those calculated with the virial theorem: 
these two values often differ by one order of magnitude, even considering different internal 
density distributions. Interestingly, the same result is found for the millimeter dust core 
distribution by MA06, who indicate two possible confining mechanisms related either to 
external pressure or to toroidal magnetic fields, excluding the first one. The hypothesis 
that these clumps, or part of them, are gravitationally bound cannot be excluded a priori if 
they are collapsing, as could be suggested by the presence of associated IRAS sources (see 
Table 3). However a particular trend in the mass ratio cannot be recognized in these cases, 
so that it remains difficult a definite conclusion about the collapsing or dispersing status of 
these clumps. 

In applying the CF algorithm to the 12 CO(1-0) line data cube, the large number of 
blended components, as well as the presence of a diffuse component in a large fraction 
of the map, make the decomposition much more difficult and the definition of "clump" 
more uncertain. Other difficulties are related to the optical depth of this line that can 
be affected by self-absorption. However, the 12 CO(1-0) observations are characterized by a 
better sampling with respect to 13 CO(2-l), approximately 1 beam and 2 beams, respectively. 
In this situation it is actually difficult to discern which line produces the most reliable slope 
for the mass spectrum. 

Even for the 12 CO (1-0) line the previous merging procedure was applied, obtaining 168 
objects from the original 275 ones. The cataloged emission in clumps accounts for the 84% 
of the total emission in the 12 CO(1-0) map. 

It is noteworthy that the positions of the 13 CO (2-1) clumps and the MA06 dust cores 
coincide with those of 12 CO (1-0) clumps, suggesting that this line, although optically thick, 
can be anyway used to reliably track the gas clumps. 

The method used to determine the 12 CO(1-0) clump masses is based on Equation 2: 
even if this empirical relationship can be questioned in the case of small-scale structures, 
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it is the only way to exploit the 12 CO(1-0) data for mass determinations and is often used 
to derive the clump mass spectra down to the lower masses (see, e.g., Heithausen et al. 
1998). The resulting mass spectrum presented in Figure 14 shows, in the case of the merged 
sample, a "flat" central part and an approximately linear decrease in the large mass limit of 
the distribution. The minimum mass and the completeness limits are, respectively, 2.9 M Q 
and 12.6 M . The spectral indexes obtained with and without merging are significantly 
different, a, — 1.3 ± 0.1 and a — 1.8 ± 0.1, respectively. While the second value is consistent 
with the corresponding index derived for 13 CO(2-l), those obtained after merging are clearly 
different. It is plausible that the high optical depth of 12 CO(l-0)makes much more difficult 
to define the actual borders of clumps. 

Note that the spectral indexes found in VMR-D for the corrected clump lists agree with 
those reported in the literature for molecular line observations (1.3 < a < 1.9, see, e.g., 
Mac Low & Klessen 2004), and are smaller than the typical slopes of the stellar IMF (e.g., 
a ~ 2.35 in Salpeter 1955). 

The slopes obtained here seem also comparable with the power-law index a^ust — 1-45 ± 
0.2 reported in MA06 if we consider the quoted errors and the fact that the dust continuum 
observations completely include those high-mass clumps which are only partially contained 
in our map (in particular the massive clump hosting IRS 16). 

The difference found between the 13 CO(2-l) and the dust indexes could be due to opacity 
effects that in the case of the more massive clumps can be considerable. This corresponds 
to an underestimate of the larger masses and then to a steeper slope in the mass spectrum. 
For another similar case in which a gas > adust, see Bains et al. (2006). 

Considering now the mass-radius relation shown in Figure 15 we can further extend 
the discussion in MA06 about the slope of this relation. In that work the exponent of the 
M oc R x relation has been evaluated x = 1.7, with this value decreasing when the masses are 
determined by adopting two different temperature-opacity values for cores with and without 
associated IRAS sources. Here on the other hand, with reference to a weak linear behavior 
in the bi-logarithmic plot, we find significantly different slopes in the case of 12 CO(1-0) and 
13 CO(2-l) clumps: x = 2.5 ±0.3 and x = 1.9 ±0.5, respectively. This discrepancy can result 
from all the uncertainties involved in the two derived mass distributions, in particular from 
the different optical depths. The 13 CO(2-l) value appears more consistent with that of MA06, 
both suggesting the presence of Bonnor-Ebert clumps. Conversely, the slope derived from 
12 CO(1-0) is more consistent with a turbulent fragmentation scenario (see, e.g., Elmegreen 
& Falgarone 1996). 

Another relevant information resulting from the cloud decomposition procedure is the 
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internal velocity dispersion a of each clump that can be used to examine its scaling relation 
with the radius R. This is generally reported as a oc R 13 (empirically derived by Larson 
1981), with (3 varying between 0.2 and 0.5 (see, e.g. Schneider & Brooks 2004), the most 
quoted value being f3 = 0.4 ±0.1 (Mac Low & Klessen 2004, and references therein). In 
Figure 16 we present the scatter plot of the velocity dispersion vs radius for the clumps 
observed in the two lines. If we fit the data with a linear regression, the resulting slope is 
(3 = 0.4 ± 0.1 for 12 CO(1-0), and (3 = 0.4 ± 0.2 for 13 CO(2-l). These values are not far 
from j3 = 0.5, that is the value corresponding to the virial equilibrium, a condition implying 
equipartition between self-gravity and turbulent energies. The same equilibrium is consistent 
with the slope x = 1.9 we find for the mass vs radius relation for the 13 CO(2-l). 

On the other hand, other authors (see, e.g., Mac Low & Klessen 2004, and references 
therein) consider a non-static scenario and explain this value with the presence of supersonic 
turbulent cascades that appear to be more consistent with the results we obtained above on 
the gravitational stability of the 13 CO(2-l) clumps. Finally, it is noteworthy that (3 = 0.5 
is also the value determined in two surveys of the galactic plane carried out by Dame et al. 
(1986) and Solomon et al. (1987). 

5. Summary 

In this work we examined the CO distribution in one of the most interesting parts of 
the Vela Molecular Ridge (i.e., a ~ 1° x 1° area of Cloud D) with unprecedented spatial 
resolution (~ 50"). The main reasons of interest for this cloud are its star formation activity 
and its location in the galactic plane. VMR-D is in fact the nearest massive star forming 
region with this galactic position. 

The millimeter line observations of the CO isotopes in VMR-D, carried out at ESO- 
SEST, revealed a complex spatial and spectral behavior down to the smallest spatial scales 
investigated (~ 1 pc). The two rotational transitions analyzed, 12 CO(1-0) and 13 CO(2-l), 
due to their different optical depth, allowed us to probe both low and high density regions, 
confirming a tight link between the gas distribution and the location of the brightest point- 
like far-IR sources known in the field. The main filamentary structures, appearing in the 
integrated intensity maps, show an arc-like shape although it is not clear if this is produced 
by expanding shells or stellar winds, both driven by nearby young (massive) stars. If an 
expanding Hn region is considered, its dynamical age is compatible with that estimated 
for the star formation activity. Contributions from nearby known supernova remnants are 
excluded. Another possibility is that the observed arcs are the product of the internal cloud 
turbulent motions. 
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The analysis of the velocity field has shown several well-separated zones, emitting at 
different velocities and suggesting a complex structure and dynamics. In two regions, the 
presence of both two velocity components (separated by ~ 3.5 and ~ 4.5 km s _1 , respectively) 
and arc-shaped morphology in the vel-pos diagrams can be interpreted as a signature of 
expanding shells. These, compressing the g be thought as the cause of the observed 

cluster formation. The evolutionary times involved in such a scenario are consistent with the 
star formation age of this cloud. Further arc-like structures are also present in the vel-pos 
diagrams which can be associated with the molecular arcs in the integrated intensity maps. 

The line profiles near the bright far-IR sources often appear clearly broadened, suggest- 
ing the presence of outflows. A search for such outflows in 12 CO(1-0) line near the IRAS 
red sources shows that, within the limits of resolution and noise of our observations, there 
are 13 candidate objects worth of further and more detailed observations for mapping their 
structure and accurately determining their physical parameters. 

One of the immediate developments of this work will be a multi-wavelength analysis of 
this cloud with the aim of correlating the clustered star formation (Lorenzetti et al. 1993; 
Massi et al. 2000, 2006a) with the gas and dust physical conditions. 

Cloud decomposition has been carried out by means of the CLUMPFIND code and 
subsequently corrected for spurious effects due to the misinterpretation of the random noise 
fluctuations in the spectra. We find that the clump mass spectra of 12 CO(1-0) and 13 CO(2- 
1) observations are different. Both values are typical for interstellar molecular clouds and, 
significantly, are in agreement with that obtained by MA06 for dust cores in VMR-D. 

In all the investigated 13 CO(2-l) clumps the masses have been found to be lower than the 
corresponding virial masses, revealing that these clumps are either gravitationally unbound 
or collapsing. 

A power-law fit of the clump mass vs radius shows that the exponent are x = 2.5 for 
12 CO(1-0) and x = 1.9 for 13 CO(2-l). Finally, a similar fit of the internal velocity dispersion 
vs clump radius gives a slope (3 ~ 0.4 in both transitions, a value suggesting either virial 
equilibrium or internal turbulence. The first hypothesis is however excluded when we take 
into account the inconsistency between observed and virial clump masses. 
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Table 1. Velocity ranges of the main emitting regions in VMR-D 



Region 


12 CO(1-0) 


13 CO(2-l) 




(km : 




(km 




North-East 


-2-: 


- 12 


OH- 4, 


6H- 12 


North-West 


0-=- 6, 10-r 18 


4H 


- 10 


Center 


4-=- 


16 


4H 


- 12 


West 





10 


0- 


h6 


South (diffuse component) 


OH 


-6 






South-East 


4H- 


16 


6H 


- 12 



Table 2. IRAS PSC sources in the observed field, with fluxes increasing from 12 /im to 

60 /im, and associated outflows. 



IRAS 


Fl2 


F 25 


*60 


^100 


LFIR Vl sr 


Outflow Wing Vrangc 


Offset 


name 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(L Q ) (kms- 1 ) 


(Y/N/. ..) (kms- 1 ) 


(") 



1 


08438-4340 (IRS 16) 


13.4 


56 


638 


1580 


566 


7.68 


Y 


bluet 
red 


6.48 
3.84 


(0, -50) 
(+50, +50) 


2 


08447-4309 


< 0.25 


< 0.4 


< 4.4 


36 


6 


11.28 










3 


08448-4343 fTTiS 17") 


8.7 


88 


327 


1010 


415 


3.96 


Y 


blue 
redt 


8.64 
8.88 


CO 0) 
(-50,0) 


4 


08448-4341 


1.3 


6.7 


< 327 


< 1010 


258 


4.20 


N 








5 


08453-4335 


< 0.34 


0.59 


< 6.4 


< 49 


9 


6.12 


N 








6 


08454-4307 


< 0.25 


0.36 


< 3.3 


< 39 


6 


2.52 


Y* 


blue 
red 


1.92 
6.72 


(+100,+50) 
(+50,0) 


7 


08458-4332 


1.1 


2.7 


18 


53 


20 


12.00 










8 


08459-4338 


< 0.29 


0.57 


6.5 


38 


8 


6.00 


Y* 


bluet 


4.32 


(-100,0) 


9 


08461-4314 


0.69 


0.79 


5.9 


< 49 


9 


2.52 










10 


08463-4343 


< 0.29 


0.49 


7.4 


< 44 


8 


10.44 


Y* 


blue 
red 


3.48 
6.24 


(+50,0) 
(+50.+100) 


11 


08465-4320 


< 0.25 


< 1.1 


< 3.5 


30 


6 


12.00 










12 


08468-4345 


0.39 


< 0.87 


< 4.8 


< 36 


7 


9.60 


Y* 


redt 


6.12 


(+100,0) 


13 


08468-4330 


< 0.3 


0.43 


3.4 


< 32 


5 


6.72 


N 








14 


08470-4321 (IRS 19) 


45 


130 


343 


407 


514 


12.96 


Y 


bluet 
red 


4.56 
6.00 


(0,-50) 
(0,-50) 


15 


08471-4346 


< 0.28 


< 0.5 


5.4 


36 


7 


8.16 










16 


08472-4326A 


0.85 


0.94 


11 


< 407 


45 


12.48 


N 








17 


08474-4325 


< 0.3 


0.96 


< 16 


58 


15 


13.20 


Y* 


bluet 
red 


8.16 
4.68 


(0,0) 
(0,+50) 


18 


08475-4352 


< 0.44 


< 0.45 


< 2.9 


31 


5 


8.52 










19 


08476-4306 (IRS 20) 


5.7 


44 


216 


504 


234 


2.28 


Y* 


blue 
red 


2.64 
4.68 


(0,0) 
(0,-50) 


20 


08477-4359 (IRS 21) 


9 


26 


317 


581 


265 


7.80 


Y* 


bluet 
redt 


5.88 
4.44 


(+50,-50) 
(0,-50) 


21 


08478-4403 


0.46 


0.51 


5.4 


< 581 


56 


10.20 










22 


08478-4303 


< 0.39 


1.6 


< 216 


< 35 


108 


4.92 


N 








2.3 


08479-4311 


< 0.32 


0.4 


< 8.7 


< 38 


8 


2.52 


Y* 


blue 
red 


3.60 
2.52 


(0,+50) 
(+50,+50) 


24 


08483-4305 


1.5 


2.5 


< 39 


189 


42 


3.36 


Y* 


blue 


7.80 


(+100,0) 



to 



Table 2 — Continued 



IRAS 


Fl2 


F 25 


F 60 


^100 


Lfir 


Vl sr Outflow 


Wing 


Vrangc 


Offset 


name 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Lq) 


(kms- 1 ) (Y/N/. ..) 




(km s- 1 ) 


(") 
















red 


2.04 


(+50,+50) 


25 08484-4350 


0.33 


< 0.56 


< 3.7 


23 


5 


6.24 Y* 


blue 


1.80 


(-50,0) 
















red 


3.24 


(-100,0) 



*A second component present in the spectrum has been preliminarily fitted and subtracted. 
tThe wing profile is probably contaminated by the contribution of very close velocity components. 



Table 3. List of the clumps detected in the 13 CO(2-l) line. 



a (2000) <5 (2000) V lsr AV Radius Mass M vir Ty,MAX J t v dV Map limit Dust IRAS 

h m s ° ' " (km s" 1 ) (km s" 1 ) (pc) (M ) (Mg) (km s" 1 )" 1 clump(s) source(s) a 



VMRD1 


8 45 


35 


-43 


52 


02 


5.9 


2.55 


0.6 


83 


842 


0.90 


2.58 


X 


MMS1,2,3 


1 


VMRD2 


8 45 


36 


-43 


48 


42 


5.3 


1.80 


0.4 


16 


248 


0.62 


0.66 


X 


MMS1,2,3 


1 


VMRD3 


8 46 


17 


-43 


56 


12 


4.1 


1.14 


0.5 


17 


142 


0.26 


0.34 








VMRD4 


8 46 


31 


-43 


54 


32 


2.9 


1.09 


0.7 


48 


165 


0.69 


0.75 




MMS4 


3,4 


VMRD5 


8 46 


35 


-43 


54 


32 


5.0 


1.53 


0.7 


70 


324 


0.97 


2.39 




MMS4,5,6 


3,4 


VMRD6 


8 46 


49 


-43 


50 


22 


4.7 


1.06 


0.4 


9 


94 


0.26 


0.06 








VMRD7 


8 46 


54 


-43 


54 


32 


4.1 


1.71 


0.6 


38 


388 


0.33 


0.17 




MMS5,6 




VMRD8 


8 47 


26 


-43 


48 


42 


5.6 


1.11 


0.4 


10 


99 


0.32 


0.46 








VMRD9 


8 47 


26 


-43 


52 


52 


2.9 


1.36 


0.4 


17 


145 


0.45 


0.44 








VMRD10 


8 47 


31 


-43 


52 


52 


4.1 


0.93 


0.4 


8 


67 


0.41 


0.39 








VMRD11 


8 47 


40 


-43 


43 


42 


11.3 


0.99 


0.3 


11 


72 


0.59 


0.50 






7 


VMRD12 


8 47 


45 


-43 


47 


02 


6.8 


0.84 


0.5 


15 


78 


0.28 


0.22 








VMRD13 


8 47 


45 


-43 


50 


22 


5.0 


1.31 


0.4 


15 


143 


0.42 


0.71 




... 

to 


8 


VMRD14 


8 47 


50 


-43 


47 


52 


5.9 


0.77 


0.4 


8 


46 


0.22 


0.21 




8 


VMRD15 


8 47 


54 


-43 


27 


52 


2.0 


0.73 


0.2 


3 


16 


0.34 


0.42 


X 


... 




VMRD16 


8 47 


59 


-43 


39 


32 


10.4 


0.64 


0.2 


2 


17 


0.22 


0.12 




MMS7 




VMRD17 


8 47 


59 


-43 


25 


22 


2.3 


0.68 


0.2 


3 


19 


0.62 


0.20 


X 




9 


VMRD18 


8 47 


59 


-43 


20 


22 


1.1 


0.50 


0.2 


3 


12 


0.78 


0.60 


X 






VMRD19 


8 48 


03 


-43 


51 


12 


4.7 


0.76 


0.2 


3 


24 


2.58 


1.44 








VMRD20 


8 48 


13 


-43 


47 


02 


7.7 


1.62 


0.3 


3 


143 


0.28 


0.16 








VMRD21 


8 48 


13 


-43 


22 


52 


2.3 


1.07 


0.4 


15 


86 


1.03 


0.73 








VMRD22 


8 48 


17 


-43 


17 


52 


1.4 


1.08 


0.2 


5 


36 


0.62 


0.62 








VMRD23 


8 48 


17 


-43 


26 


12 


2.3 


0.88 


0.2 


3 


23 


0.57 


0.65 








VMRD24 


8 48 


27 


-43 


41 


12 


6.5 


1.69 


0.3 


5 


201 


0.48 


0.11 








VMRD25 


8 48 


31 


-43 


22 


02 


2.0 


0.85 


0.2 


4 


26 


0.54 


0.44 








VMRD26 


8 48 


31 


-43 


25 


22 


2.3 


0.91 


0.2 


6 


34 


0.42 


0.53 








VMRD27 


8 48 


31 


-43 


34 


32 


12.5 


1.18 


0.4 


11 


125 


0.47 


0.24 






11 


VMRD28 


8 48 


40 


-43 


29 


32 


11.3 


1.33 


0.4 


8 


153 


0.50 


0.45 




MMS8,9,12,14,15,16 


14 


VMRD29 


8 48 


40 


-43 


17 


52 


2.3 


0.87 


0.2 


6 


23 


0.88 


0.53 








VMRD30 


8 48 


50 


-43 


29 


32 


13.1 


1.07 


0.5 


16 


123 


0.39 


0.20 




MMS8,9,12,13,14,15,16 


14 


VMRD31 


8 48 


54 


-43 


15 


22 


2.0 


0.97 


0.2 


8 


35 


0.79 


0.78 








VMRD32 


8 48 


59 


-43 


59 


32 


8.3 


1.44 


0.5 


16 


229 


0.81 


0.53 






15 


VMRD33 


8 49 


03 


-43 


32 


52 


11.9 


0.86 


0.2 


4 


37 


0.58 


0.54 









Table 3 — Continued 





ot (ZUUUJ 


X /onnn'i 
{ZVDV) 


VI sr 


A \f 
L\ V 


Radius 


Mass 




"TV, MAX 


J tv dV 


Map limit 


Dust 






h m s 


o / // 


(Km s ) 


(Km s ) 


(P c ) 


( M 0j 


l M 0) 


(Km s ) 






clump(s) 


source 


VMRD34 


8 49 03 


-43 37 52 


11.9 


1.69 


0.2 


10 


143 


0.54 


0.71 




MMS17,18 


16 


VMRD35 


8 49 12 


-43 22 02 


1.4 


1.27 


0.3 


5 


93 


0.37 


0.52 








VMRD36 


8 49 13 


-43 36 12 


11.3 


1.55 


0.3 


32 


170 


0.94 


1.81 




MMS19,20,21 


17 


VMRD37 


8 49 13 


-43 36 12 


12.8 


0.86 


0.3 


9 


45 


0.64 


0.46 




MMS19,20,21 


17 


VMRD38 


8 49 22 


-44 00 22 


8.3 


1.40 


0.6 


19 


228 


0.39 


0.24 








VMRD39 


8 49 22 


-43 58 42 


7.4 


1.07 


0.4 


8 


88 


1.01 


0.81 








VMRD40 


8 49 26 


-43 22 02 


3.5 


0.72 


0.2 


3 


21 


0.32 


0.30 








VMRD41 


8 49 26 


-43 17 02 


2.6 


0.90 


0.3 


10 


59 


0.28 


0.37 




MMS22,24 


19 


VMRD42 


8 49 26 


-43 17 02 


1.7 


0.85 


0.2 


3 


26 


0.26 


0.23 




MMS22,24 


19 


VMRD43 


8 49 27 


-44 10 22 


8.9 


1.57 


0.7 


78 


365 


3.55 


5.50 




MMS25, 26,27 


20 


VMRD44 


8 49 31 


-43 20 22 


7.7 


1.16 


0.3 


4 


73 


0.40 


0.42 








VMRD45 


8 49 31 


-44 07 52 


9.8 


1.63 


0.7 


141 


373 


1.10 


0.59 




MMS23.25, 26,27 


20 


VMRD46 


8 49 49 


-43 21 12 


8.3 


1.01 


0.3 


9 


66 


0.63 


0.69 






23 


VMRD47 


8 49 59 


-44 00 22 


6.8 


1.00 


0.2 


2 


36 


0.32 


0.30 






. .to 

oo 


VMRD48 


8 50 03 


-43 22 52 


9.5 


1.27 


0.4 


20 


136 


0.60 


0.81 


X 




VMRD49 


8 50 07 


-43 16 12 


4.1 


0.88 


0.2 


4 


28 


0.35 


0.30 


X 


MMS28 


24"' 



a The identification numbers are the same as in Table 2. 




Fig. 1. — Typical line profiles. This sample illustrates the part of the map corresponding to 
the offset interval +750" < a off < +1150", +1350" < S off < +1700". 
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5 10 15 20 
V (km s" 1 ) 



Fig. 2. — Sum of all the 12 CO(1-0) (upper panel) and 13 CO(2-l) spectra (lower panel). The 
dotted line in both panels represents the contribution of the sole volume pixels assigned 
to clumps after the cloud decomposition (see Section 4). Since the solid lines can be also 
interpreted as averaged spectra (by dividing by the number of points showing significant 
emission), a re-scaled y axis is reported on the right to help in this reading. 




1000 500 -500 -1000 -1500 -2000 



Fig. 3. — Integrated intensity map for the 12 CO(1-0) emission in VMR-D region, in the range 
—2 20 km s _1 . Contour levels start from 5 K km s -1 and are separated by 15 K km s _1 . 
Dashed lines delimit the observed area. On the right, the spatial scale corresponding to the 
estimated distance (d = 700 pc) is shown. 
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Fig. 5. — Superposition of the 13 CO(2-l) contour levels (shown in Figure 4) on the 12 CO(l- 
0) integrated intensity map (shown in Figure 3). The positions of the "red" (F±2 < F25 < Fqq) 
IRAS sources are also marked with filled circles. 
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Fig. 6.— Channel maps from -2 to 8 km s" 1 for the 12 CO(1-0) (left panels) and 13 CO(2- 
1) (right panels) emission in VMR-D region, in the velocity ranges indicated on the left 
border of each box. For 12 CO(1-0) Contour levels start from 2 K km s _1 and are sepa- 
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5 10 1 5 20 5 10 15 20 5 10 15 20 5 10 15 20 5 10 15 20 5 10 15 20 

VJkm s 1 ) 

Fig. 8. — Vel-Soff diagrams of the 12 CO(1-0) emission in VMR-D, for fixed values of cc //, 
indicated in each panel. Contour levels start from T m ^ = 2 K and are separated by 4 K. Four 
arc-like structures have been recognized and marked with a solid arc of ellipse, each one in 
the diagram in which it is best visible. 
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Fig. 9. — A 2° x 2° field centered on the (0,0) position of the map is shown, with the contours 
of the 13 CO(2-l) emission plotted as in Figure 4. The bold dotted circles roughly indicate 
three arc-like emitting regions, while the thin dashed lines mark the border of the observed 
zone. Crosses indicate the center of the Hn regions; for two of them with significant spatial 
extent, Gum 18 and RCW 35, an ellipse is also drawn whose axes length is taken from 
literature. The remaining Hll regions are compact objects and cannot be responsible for 
the large-scale morphology. The locations of OB-type stars (filled circles), with their proper 
motion track extrapolated for the last 10 5 yr, are also marked. 
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Fig. 10. — Plot of the dynamical age td yn (solid line for d = 700 pc, and dotted lines for 
d = 500 pc and d = 900 pc, respectively) for a Hll region and of the time at which the 
fragmentation of the driven shell starts (for different values of the sound velocity c s and a 
radius of ~ 3.4 pc) as a function of the density n . 
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Fig. 11. — Panel a: 12 CO(1-0) spectrum corresponding to the closest pointing to the IRS 17 
location, at map offset {a fj = — 1350", r) // — 500"); the dashed line is the gaussian deter- 
mined as in WB99. Panel b and c: interpolated greyscale maps for blue and red outflow 
components, respectively. Greyscale and contour levels are the same in both maps; the levels 
are in steps of 5 K km s -1 . Observed positions (crosses) and IRS 17 location (filled circle) 
are also shown. The blue wing peak assigned to the outflow is coincident with the central 
position in panel b, while the brighter feature at 100" on the right is considered unrelated. 
Note that the red wing is probably contaminated by the contribution of a second component 
to the line profile. Panels d, e, /: the same as panels a, b, c, but for IRAS 08474-4325 (even 
here the blue wing is probably contaminated). 
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Fig. 12. — Example of bad clump detection performed by CF. Panel a: a grid of 12 CO(l- 
0) spectra (resampled in resolution from 0.12 to 0.6 km s" 1 ) around two positions (marked 
with a thicker line) recognized by CF as centroids of different clumps, with almost coincident 
peak velocity (Vi sr ~ 7 km s" 1 ). Because these components satisfy our merging criteria 
(see text), we conservatively chose to consider them as a single clump. Note that a second 
component at V/ sr ~ 10 km s -1 is clearly resolved in some spectra and is correctly assigned by 
CF to another clump. Panel b: the two thicker spectra shown in panel a are superimposed 
and plotted with solid and dashed line, respectively. The velocity channels erroneously 
assigned by CF to two different clumps are highlighted (bold line). 
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Fig. 13. — Clump mass spectrum of VMR-D derived from 13 CO(2-l) observations using 
CF. The dashed vertical lines indicate the minimum mass (leftmost) and the completeness 
(rightmost) limits, respectively. The dotted line shows the CF mass spectrum while the solid 
line is obtained after the merging procedure (see text). The linear best-fit is also shown in 
both cases, whose slope corresponds to a — 1, where a is the exponent of the power law 
d.Y.'d.U oc M n . 



-42- 




Fig. 14.— As in Figure 13, but for 12 CO(1-0). 
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Fig. 15. — Bi-logarithmic plot of the mass vs radius relation for the clumps resulting from the 
supervised clump decomposition of VMR-D, for the 12 CO(1-0) and the 13 CO(2-l) transitions, 
respectively. In both panels, the dashed line indicates the linear trend of the data. The slopes 
are x = 2.5 ± 0.3 for 12 CO(1-0), and x = 1.9 ± 0.5 for 13 CO(2-l). 
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Fig. 16. — Bi-logarithmic plot of the internal velocity dispersion vs radius for the VMR-D 
clumps. The data are arranged as in Figure 15. The slopes are f3 = 0.4 ±0.1 for 12 CO(1-0), 
and (3 = 0.4 ± 0.2 for 13 CO(2-l). 



